Lateral correlation of multivalent counterions is the universal mechanism of charge 

inversion 
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We review works on screening of a macroion, such as colloidal particle or double helix DNA, by multivalent 
counterions. Multivalent metal ions, charged micelles, short or long polyelectrolytes can play the role of multivalent 
counterions. Due to the strong Coulomb repulsion such multivalent counterions form a strongly correlated liquid 
resembling a Wigner crystal at the surface of the macroion. Even if the macroion is neutralized by this liquid, a 
newly arriving counterion creates in the liquid a correlation hole or image which attracts the ion to the surface. 
As a result, total charge of adsorbed counterions can substantially exceed the bare macroion charge, so that the 
macroion net charge inverts its sign. We discuss two previously suggested explanations of charge inversion and 
show that physics underlying both of them is intrinsically that of correlations, so that correlation is the universal 
force driving charge inversion. 



I. INTRODUCTION 



charged surfaces was recognized even earlier 



o 

CO 
C3 



T3 
£ 
O 

> 

m 
O 



O ' 

o : 

-i— > ■ 

a ■ 



i 

-a 
c 

o 
o 



X 



Charge inversion is a counterintuitive phenomenon in 
which a strongly charged particle (a macroion) binds so 
many counterions that its net charge changes sign. As 
shown below the binding energy of a counterion with 
large charge Z is larger than UbT, so that this net charge 
is easily observable: it is the net charge that determines 
a particle drift in a weak field electrophoresis. Charge 
inversion is possible for a variety of macroions, ranging 
from the charged surface of mica to charged lipid mem- 
branes, colloids, DNA or actin. Multivalent metal ions, 
small colloidal particles, charged micelles, short or long 
polyelectrolytes including DNA can play the role of mul- 
tivalent counterions. Recently charge inversion has at- 
tracted significant attention 

Theoretically, charge inversion can be also thought of 
as an over-screening. Indeed, the simplest screening at- 
mosphere, familiar from the linear Debye-Hiickel (DH) 
theory, compensates at any finite distance only a part 
of the macroion charge. It can be proven that this is 
true also in the non-linear Poisson-Boltzmann (PB) the- 
ory. The statement that the net charge preserves sign of 
the bare charge agrees with the common sense. One can 
think that this statement is even more universal than 
results of PB equation. However, this presumption of 
common sense fails for screening by Z-valent counterions 
with large Z (Z-ions). In this case, most of counteri- 
ons are localized at the very surface of macroion. The 
energy of their lateral Coulomb interaction may exceed 
kgT by an order of magnitude or more. As a result HQ, 
at the macroion surface ZAons form a two-dimensional 
(2D) strongly correlated liquid (SCL) with properties re- 
sembling a Wigner crystal (WC), which is shown in Fig. 
1. 

The negative chemical potential of this liquid leads to 
a correlation induced attraction of Z-ions to the surface, 
which brings more of them to the surface than neces- 
sary to neutralize it. Role of correlations of Z-ions in 
another physical phenomenon attraction of two likely- 
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FIG. 1. The Wigner crystal of positive Z-ions on the neg- 
ative uniform background of surface charge. A hexagonal 
Wigner-Seitz cell and its simplified version as a disk with ra- 
dius R are shown. 

This paper consists of two parts. In the first part, we 
review recent works on correlation induced charge inver- 
sion. Sec. II deals with small Z-ions, while in Sec. Ill 
the role of Z-ions is played by polyelectrolytes. In the 
second part (Sec. IV), we discuss proposed alternative 
mechanisms of charge inversion: metallization approach 
(mental smearing of Z-ions on the macroion surface) and 
counterion release. It is shown there that they are also 
based on the correlation physics, so that some of the ap- 
parent difference is purely semantic. 



II. CORRELATIONS OF SMALL SIZE 
MULTIVALENT COUNTERIONS 



Let us consider screening a macroion surface with a 
negative charge density —a by Z : 1 salt with positive 
Z-ions in the case of large enough Z and a. If such 
a surface is neutralized by Z-ions, the average distance 
between them in directions parallel to the plane equals 
Rq = {-kct / Ze) -1 / 2 . If this distance is much larger than 
the Gouy-Chapman length A = DksT '/ '(2-KZea), which 
is the typical distance to the surface mean field PB ap- 
proach obviously fail. The ratio of these two lengths 
Ra/X = 2r, where V = Z 2 e 2 /(R Dk B T) is the dimen- 
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sionless inverse temperature measured in units of the typ- 
ical interaction energy and D is the dielectric constant of 
water. Here we are concerned with large enough Z and 
cr at which parameter r> 1. For example, at Z = 3 and 
a = 1.0 e/nm 2 we get T = 6.4, A ~ 0.08 nm and R ~ 1.0 
ran. In such conditions Gouy-Chapman solution of PB 
equation can not be valid any more. Indeed, if distance 
of a Z-ion to the surface x is in the interval A -C x -C Rq, 
it does not feel other Z-ions but interacts only with the 
macroion surface. This interaction leads to an exponen- 
tial decay of concentration of Z-ions as function of x in 
the range x <C Ro instead of 1/x 2 Gouy-Chapman law. 
A new theory based on inequality r>l was suggested in 
Ref . [HJJ] • The main idea of this theory is that at T > 1 
the screening atmosphere can be divided in two distinct 
phases: 2D SCL, which should be treated exactly, and 
the gas phase which at large distances x can be desribed 
by the PB equation and is in equilibrium with 2D SCL. 

Thermodynamic properties of 2D liquid of classical 
charged particles on the neutralizing background, so 
called one-component plasma (OCP) are well known (see 
bibliography in Ref. ||]4|). The chemical potential of 
OCP can be written as /i = fj,id + /i c , where [lid and fi c 
are its ideal and correlation parts. At T — OCP forms 
WC. To calculate fi c , we start from the energy of WC 
per Z-ion, 

e(n) ~ -l.llZ 2 e 2 /RD = -1.96n 1/2 Z 2 e 2 /D, (1) 

where R — (7m) _1 / 2 is the radius of its Wigner-Seitz 
(WS) cell (See Fig. 1) and n is 2D concentration of Z- 
ions. This gives // c = ^wc, where 

mc = feM = _i. 6 5rfc B T = -\MZ 2 e 2 /DR. (2) 
on 

In the range of temperatures where 1 -C V -C 130 WC is 
melted and for SCL \jl c = (Awe + where 8^ is a small 
positive correction H,||. 

Physics of a large and negative Hwc can be understood 
as follows. Let us imagine for a moment that an insulat- 
ing macroion is replaced by a neutral metallic particle. In 
this case, each Z-ion creates an image charge of opposite 
sign inside the metal. Energy of attraction to the image 
is U(x) = —{Ze) 2 /ADx, where x is the distance from Z- 
ion to the surface. At the metal surface energy of the 
interaction with the image equals U(a) = —(Ze) 2 /4:Da, 
where a is radius of the countcrion. The chemical poten- 
tial fiwc f° r a charged surface of an insulating macroion 
can be interpreted in a similar language of images. 



FIG. 2. The origin of attraction of a new positive Z-ion 
to the already neutralized surface. Z-ions are shown by solid 
circles. The new Z-ion creates its negative correlation hole. 

Consider bringing a new Z-ion to the macroion surface 
already covered by an adsorbed layer of Z-ions (See Fig. 
2). This layer plays the role of a metal surface. Indeed, 
the new Z-ion repels nearest adsorbed ones, creating a 
correlation hole for himself. In other words, it creates 
a negative image. Calculation of the energy of attrac- 
tion to the image in this case is, however, less trivial 
than in the case of a metal. The problem is that min- 
imal size of the image in the adsorbed layer is equal to 
the WS cell radius R. (The adsorbed layer is a good 
metal only at larger scales.) Thus, for WC we arrive at 
Hwc ~ —{Ze) 2 /DR. Eq. (g) provides the numerical co- 
efficient in this expression. It is clear now that contrary 
to the case of a metallic surface, charge inversion for in- 
sulating macroion requires a finite density of an adsorbed 
layer, and, correspondingly, a finite bare charge density 
cr. 

Correlation induced attraction of Z-ion to a neutral 
SCL has two interesting analogies in the solid state and 
atomic physics. The energy |/LtvKc| plays the same role 
for Z-ions as the work function plays for electrons of a 
metal or the ionization energy plays in a many-electron 
atom. It is known that in the Thomas-Fermi approx- 
imation which neglects exchange and correlation holes, 
both the work function of a meta l PjJ and the ionization 
energy of a many-electron atom |25| vanish. For screen- 
ing by Z-ions the PB approximation is an analog of the 
Thomas- Fermi approximation and it results in /i c = 0. 
Only correlations lead to a finite work function or ioniza- 
tion energy and to a finite \i c ~ Hwc- 

Let us now discuss behavior of the concentration of 
Z-ions, N(x) near the surface. For this purpose let us 
extract a Z-ion from SCL and move it along the x axis. 
As it is mentioned above, this Z-ion leaves behind its 
correlation hole. In the range of distances x <C R, the 
correlation hole is a disc of the surface charge with radius 
R (WS cell) and the Z-ion is attracted to the surface by 
its uniform electric field E = 2tt<t/D. We can say that 
our Z-ion does not see other Z-ions until x <C R. There- 
fore, N(x) = N s exp{-x/\) at x < R. Here N s = n/X 
is the three-dimensional concentration of Z-ions at the 
plane. At large enough x the concentration N(x) should 
saturate at the level of 

N Q = N s exp(-\»wc\/k B T) = iV s exp(-1.65r). (3) 



Let us find out how and where this saturation happens. 
At x ~ R = 2TX the ratio x/X reaches roughly speak- 
ing a half of \nwc\/kBT. At x ^> R the Z-ion cre- 
ates the image with relatively small radius R which at- 
tracts Z-ion back and provides the decaying Coulomb 
correction to the activation energy of N(x): N(x) = 
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N s exp (-[\nwc\ - Z 2 e 2 /4Dx]/k B T). This correction is 
similar to the well known " image" correction to the work 
function of a metal @. At x = Z 2 e 2 /ADk B T = \T 2 /2 
this correction reaches fc^T, so that iV(x) saturates at 
the value N . 

The dramatic difference between the described above 
exponential decay of N(x) and the Gouy-Chapman 
1/(A + a;) 2 -law is obviously related to the correlation ef- 
fects. Recently such an exponential decay of N(x) with 
the following tendency to saturation was rederived in 
more formal way and confirmed by Monte-Carlo simu- 
lations p6[| . At distances x XT 2 /2, interaction of the 
removed ion with its correlation hole in SCL is not impor- 
tant and the correlation between ions of the gas phase are 
even weaker because N(x) is exponentially small. There- 
fore, at larger distances one can describe N(x) by PB 
equation J^fl]. 

Studying charge inversion, we want to calculate the net 
charge density of the macroion surface a* = —a + Zen > 
0. Even though the system is not neutral, we can still 
use the chemical potential of a neutral OCP given by 
Eq. (|J) after the following exact transformation. Indeed, 
let us add uniform charge densities —a* and a* to the 
macroion surface. The first addition results in a neutral 
OCP. The second addition creates only uniform poten- 
tial eip(0) on the macroion surface. For example, if the 
macroion is a sphere with radius r and screening radius 
of solution is larger than r we get ip(0) — Q* /Dr, where 
Q* = 47rr 2 er* = —Q + Airr 2 neZ is the net charge of the 
sphere and —Q = 47rr 2 er is its bare charge. It is impor- 
tant to emphasize that macroscopic net charge Q* does 
not interact with OCP, because potential of the uniformly 
charged sphere at the neutral OCP is constant. 

The condition of balance of the electrochemical poten- 
tial at the surface of macroion, u, and in the bulk of 
solution, at now can be written as Zeip(0) = —^wc + 
(lib - Hid) = ~Hwc + k B T\n{N s /N) = k B Tln{N/N ), 
where N is the concentration of Z-ions in the bulk of 
solution. It is clear that when N > No, the net charge 
density a* is indeed positive, i.e. has the sign opposite 
to the bare charge density —a. The concentration No is 
very small because \nwc\/k B T = l.QbT » 1. Therefore, 
it is easy to achieve charge inversion increasing N. At 
large enough N we have — fad] <C |mvkc|- This gives 
a simple equation for calculation of the maximal inverted 
charge density 



Q* = OM^/QZe. 



(5) 



^(0) = \uwc\/Ze. 



(4) 



Let us consider a sphere with charge —Q screened by 
Z : 1 salt with a large concentration, N. In this case 
-0(0) = Q* I Dr and Eq. ([I]) has a simple meaning: 
\Hwc\/Ze is the "correlation" voltage which charges a 
spherical capacitor. Expressing R and |/zvpc| through 
Q and Z we arrive at the simple prediction for the 
maximum possible inverted charge: 



This charge is much larger than Ze, but still is smaller 
than Q because of limitations imposed by the large charg- 
ing energy of the macroscopic net charge. For example, 
for Q = lOOe, Z = 4, we get Q* = lie. Eq. (|) was 
recently confirmed by numerical simulations [T^j22 |. 

It was also shown in Ref. 0,0]) that m the case of 
a cylinder, the conventional picture of nonlinear screen- 
ing known as the Onsager-Manning (OM) condensation 
should be strongly modified when dealing with multi- 
valent ions. Consider a cylinder with a negative lin- 
ear charge density — rj and assume that rj > rj z , where 
Vz = k B TD/Ze. OM theory @, shows that such a 
strongly charged cylinder is partially screened by coun- 
terions residing at its surface, so that net linear charge 
density of the cylinder, n* = —r/ z . The rest of the charge 
is screened at much larger distances according to the lin- 
ear DH theory. The OM theory uses PB mean field ap- 
proach and, therefore, does not take into account lateral 
correlations of counterions. It is shown in Ref. |3|,[I| that 
the correlation induced negative chemical potential [iwc 
leads to inversion of the sign of rf at N > Nq in the 
case of cylinder, too. At large enough N inverted charge 
density rf , can reach k B TD/e. 

Even stronger charge inversion of a spherical or cylin- 
drical macroion can be obtained in the presence of sub- 
stantial concentration of monovalent salt, such as NaCl, 
in solution. Monovalent salt screens long range Coulomb 
interactions stronger than short range lateral correla- 
tions between adsorbed Z-ions. Therefore, screening by 
monovalent salt diminishes the charging energy of the 
macroion much stronger than the correlation energy of 
Z-ions. As a results, the inverted charge Q* becomes 
larger than that predicted by Eq. ([|) and scales linearly 
with Q. Since, in the presence of a sufficient concen- 
tration of salt, the macroion is screened at the distance 
smaller than its size, the macroion can be thought of as 
an over-screened surface, with inverted charge Q* pro- 
portional to the surface area. In this sense, overall shape 
of the macroion and its surface is irrelevant, at least to a 
first approximation. Therefore, we consider here a sim- 
ple case: screening of a planar macroion surface with a 
negative surface charge density —a by solution with con- 
centration N of Z : 1 salt and a large concentration iVi of 
a monovalent salt. Correspondingly, we assume that all 
weak interactions are screened with DH screening length 
r s = (8irl B Ni) 1 . For simplicity, we discuss here only 
a charged surface of insulating macroion with the same 
dielectric constant as water. Complications related to 
the difference between dielectric constants are discussed 
in Ref. @. 

The dependence of the charge inversion ratio, a* / 'a, on 
r s was calculated analytically [|l3| in two limiting cases 
r s Rq and r s <§C Ro, where i?o = (ira / Ze) -1 / 2 is the 
radius a WS cell at the neutral point n = cr/Ze. At 
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r s Rq calculation starts from Eq. Electrostatic 
potential can be calculated as potential of the plane with 
the charge density a* screened at the distance r s . This 
gives ip(0) = 47rer*r s . At r s 3> Rq screening by monova- 
lent ions does not change Eq. (Q) substantially so that we 
still can use it in Eq. (Q) which now describes charging of 
a plane capacitor by voltage \^wc\/Ze. This gives Q 

a* /a = 0Al(Ra/r s ) « 1 (r, » Rq). (6) 

Thus, at r s 3> Rq inverted charge density grows with de- 
creasing r s . Now we switch to the case of strong screen- 
ing, r s <C Rq. It seems that in this case a* should de- 
crease with decreasing r s , because screening reduces the 
energy of SCL and leads to its evaporation. In fact, this 
is what eventually happens. However, there is a range 
where r s <^ R , but the energy of SCL is still large. In 
this range, as r s decreases, the repulsion between Z-ions 
becomes weaker and makes it easier to pack more Z-ions 
on the plane. Therefore, a* continues to grow with de- 
creasing r s . At r s <C Rq it is convenient to minimize 
directly the free energy of WC of Z-ions at the macroin 
surface with respect of n. This free energy consists of 
nearest neighbor repulsion energies of Z-ions and the at- 
traction energy of Z-ions to the charge surface. All in- 
teractions are screened according to DH theory, so that 
interaction of non-nearest neighbors can be neglected. 
This gives || 

F = -4nar s Zen/D + (3nZ 2 e 2 /DA) exp(-A/r s ), (7) 

where A = (2/v / 3) 1 / 2 w~ 1 / 2 is the lattice constant of the 
hexagonal WC. Minimizing this free energy with respect 
to n one arrives at 

a* /a = (V2V3)[i?oAi HR /r s )] 2 (r s « Rq). (8) 

Thus a* I a grows with decreasing r s and can become 
larger than 100%. At r s ~ R Eq. (§) and Eq. (§) 
match each other. 

III. CORRELATIONS OF ADSORBED 
POLYELECTROLYTE MOLECULES 

A. Adsorption of rod-like polyelectrolytes on the 
charged plane 

A practically important class of multivalent ions are 
polyelcctrolytc (PE) molecules. In this section we discuss 
charge inversion caused by adsorption of long rod-like Z- 
ions. To make signs consistent with the case of DNA, we 
assume that PE charge is negative, — rj per unit length, 
while the macroion surface is a plane with positive charge 
density a. The basic idea is once again the strong lateral 
correlations: due to the strong lateral repulsion, charged 
rods adsorbed at the surface tend to be parallel to each 



other and have a short range order of an one-dimensional 
WC in the perpendicular to rods direction (Fig. 3). 




FIG. 3. Rod-like negative Z-ions such as double helix DNA 
are adsorbed on a positive uniformly charged plane. Strong 
Coulomb repulsion of rods leads to one-dimensional crystal- 
lization with lattice constant A. 

It is well known that the bare charge density of DNA, 
— ?7o, is about four times the critical density — r\ c = ksT/e 
of OM condensation. Therefore, about three quarters of 
the bare charge of an isolated DNA is compensated by 
positive monovalent ions residing at its surface so that 
the net charge of DNA in the bulk solution is rf = —i] c . 
This is generally not true for adsorbed DNA. Indeed, 
some OM-condensed ions, as they are repelled from the 
charged surface, can be pushed away from DNA and 
released into the solution when DNA is adsorbed ||. 
This effect is particularly important when monovalent 
salt concentration is low, otherwise counterions do not 
gain enough entropy to justify the release. 

It was shown Jl3[ that the criterion of counterion re- 
lease is governed by the comparison of screening length, 
r s , of monovalent salt and the spacing between critically 
charged rods, Ao = r] c /a, when they neutralize the plane. 
If screening is strong and r s <C Aq, the potential of the 
surface is so weak that counterions are not released. In 
other words, charge density for each adsorbed DNA he- 
lix remains r/* = —r) c . Simultaneously, at r s <C Aq the 
DH approximation can be used to describe screening of 
the charged surface by monovalent salt. Using these sim- 
plifications one can directly minimize the free energy of 
one-dimensional crystal of DNA rods on the positive sur- 
face written similarly to Eq. (Q). Then the competition 
between attraction of DNA helices to the surface and the 
repulsion of the neighboring helices results in the nega- 
tive net surface charge density — a* . The charge inversion 
ratio reads ||[l3| : 

a* I a = {ri c /(7r s )/\\i(rj c /(jr s ) {r\ c jar s > l,r s < Aq). 

(9) 

Thus the inversion ratio grows with decreasing r s as in 
the case of spherical Z-ions. At small enough r s and cr, 
the inversion ratio can reach 200% before DNA molecules 
are released from the surface. It is larger than for spher- 
ical ions, because in this case, due to the large length 
of DNA helix, the correlation energy remains large and 
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WC-like short range order is preserved at smaller val- 
ues of a and/or r s . In the works [ ^3|Jl^ ], we called this 
phenomenon "giant charge inversion." 

Let us switch now to the opposite extreme of weak 
screening by a monovalent salt, r s » Aq. In this case, 
screening of the overcharged plane by monovalent salt 
becomes strongly nonlinear, with the Gouy-Chapman 
screening length A* = DksT j (2irecr*) much smaller than 
r s . Furthermore, some counterions are released from 
DNA upon adsorbtion. As a result the absolute value 
of the net linear charge density of each adsorbed DNA, 
77*, becomes larger than rj c . Two nonlinear equations for 
a* and 77* are derived in Ref. Jl3| and are discussed below 
in Sec. IV. Their solution at r s ^> Aq reads: 



a* I a = (ri c /27racr) exp (-^ln{r s /a) ln(A /27ra) , (10) 



(11) 



77* = rj c ^J\Yi(r s /a)/ ln(A /27ra) 



At r s ~ Aq we get rj* ~ rj c , A* ~ r s and a* /a ~ 
r] c / (2irr s a) so that Eq. ( |To| ) crosses over smoothly to the 
strong screening result of Eq. Since 7/* can not be 
smaller than 77c the fact that 77* ~ 7y c already at r s ~ Ao 
proves that at r s <C Aq, indeed, rf ~ rj c . Simultaneously, 
the fact that at this point A* ~ r s means that r s <C A* at 
r s <C Aq, so that screening in this regime becomes linear 
and is described by DH theory. 



B. Flexible polyelectrolytes wrapping around the 
charged spheres 




FIG. 4. A PE molecule winding around a spherical 
macroion. Due to Coulomb repulsion, neighboring turns 
he parallel to each other. Locally, they resemble a 
one- dimensional Wigner crystal with the lattice constant A. 



The last example of the correlation driven charge in- 
version we want to mention here is the complexation in a 
solution with given concentrations of a long flexible neg- 
ative PE and positive spherical particles such as colloids, 
micelles, globular proteins or dendromers |lS[ ]. PE binds 
spheres winding around them. If the total charge of PE 
in the solution is larger than the total charge of spheres, 
repulsive correlation of PE turns on a sphere surface lead 
to inversion of the net charge of each sphere. Negative 
spheres repel each other and form on the PE molecule 
periodic beads-on-a-string structure, which resembles 10 
nm fiber of chromatin (Fig. 5). 
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Until now we talked only about adsorption of a rod- 
like PE. If the persistence length L p is finite, PE is re- 
leased from the charged surface when r s is very small. 
This happens when interaction of a PE segment of the 
length L p with the surface becomes smaller than ksT. 
On the other hand, at larger r s a flexible PE lays flat 
at the surface. Due to its electrostatic rigidity neighbor- 
ing molecules form WC-like SCL and behave similarly to 
that of a rod-like PE. The correlation induced adsorption 
of a flexible weakly charged PE on an oppositely charged 
surface was comprehensively studied in Ref. |flSfl . When 
surface is covered by less than one complete layer of PE, 
results are indeed close to presented above. 

Let us also mention the charge inversion in the problem 
of complexation of a positive sphere and a flexible neg- 
ative PE (Fig. 4). Refs. [§|§|nj predicted substantial 
charge inversion in this case: more PE is wound around 
a sphere than necessary to neutralize it. Role of corre- 
lations in such charge inversion was recently emphasized 
in Ref. |T(| . It was shown that neighbouring turns repel 
each other and form almost equidistant solenoid, which 
locally resembles WC. The tail of PE repels PE adsorbed 
at the surface and creates a correlation hole, which at- 
tracts the tail back to the surface. 



FIG. 5. Complexation of a negative PE molecule and many 
positive spheres in a necklace-like structure. On the surface 
of a sphere neighboring PE turns form WC similar to Fig. 
4. At larger scale, charged spheres repel each other and form 
one- dimensional Wigner crystal along the PE molecule. 



If the total charge of PE is smaller than the total charge 
of spheres, the latter are under-screened by PE and their 
net charges are positive. Bound spheres once more re- 
pel each other and form a periodic necklace. Because a 
segment of PE wound around a sphere interacts almost 
exclusively with this sphere, it plays the role of WS cell. 
These WC-like correlations lead to inversion of charge 
of a PE molecule: spheres bind to PE in such a great 
number that the net charge of the PE molecule becomes 
positive. It is shown in Ref. |l^] that inverted charge 
of PE by the absolute value can be larger than the bare 
charge of PE even in the limit of very weak screening by 
monovalent salt. When a large concentration of a mono- 
valent salt is added charge inversion can reach giant pro- 
portions. This theory is in a qualitative agreement with 
recent experiments on micelles-PE systems pT| . 
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IV. ARE THERE OTHER MECHANISMS OF 
CHARGE INVERSION? 

However briefly, we completed the review of works on 
correlation induced charge inversion. It is now time to 
ask: Are there alternative, correlation- independent, elec- 
trostatic mechanisms leading to this phenomenon? Our 
answer is no, and we argue that correlations-based mech- 
anism is the universal one. Specifically, we consider below 
two theories suggested in literature. One, pioneered by 
, considers adsorbed Z-ions as a smeared continuous 
medium similar to a metal. We argue that this metalliza- 
tion approach is an approximation for correlation theory; 
we discuss when this anzatz is accurate, and when it fails. 
The other theory, put forward by 0, views counterion 
release as the driving force behind charge inversion. To 
this end, we show that while counterion release is obvi- 
ously favorable for charge inversion, it is itself driven by 
correlations. 



A. Metallization approach 

This approach was clearly formulated for the cases of 
complexation of a PE with a sphere in |5| and for adsorp- 
tion of flexible polymers on a charged plane in 0. The 
physical basis of this theory is a single, seemingly be- 
nign approximation. Namely, it considers the adsorbed 
Z-ions, at the macroion surface, as a smeared continuum, 
while still treating bulk solution as consisting of discrete 
charges. Why is smearing so much favorable as to cause 
adsorbtion of Z-ions to the macroion which is already 
neutralized or even overcharged? 

The answer is that smearing amounts to complete ne- 
glect of self-energy e se if — Z 2 e 2 /2Da of adsorbed Z-ions; 
better to say, it neglects the energy of each Z-ion's elec- 
tric field in the range of distances between the ion radius 
a and macroion size r. On the other hand, correlations 
can be viewed as physical mechanism eliminating a part 
of self-energy of each Z-ion corresponding to the field in 
the range of distances between correlation hole size R 
(which plays the role of screening radius of SCL) and 
macroion size r. In this sense, metallization overesti- 
mates the role of correlations and the charge inversion 
ratio, particularly when R 3> a. 

To appreciate the difference between correlated and 
smeared set of Z-ions, it is worth comparing both to the 
random distribution of the same ions. Correlated config- 
uration is more favorable than random: correlations hap- 
pen because ions reconfigure themselves non-randomly to 
gain some free energy. Smeared continuum is also more 
favorable than random configuration of discrete ions, the 
difference being precisely the self-energy of ions. 

Let us illustrate these ideas in the simplest geometry - a 
sphere with charge — Q and radius r screened by spherical 
Z-ions of radius a <C r in the absence of monovalent salt. 



Correlation theory, as we have seen, results in Eq. (g) 
for this problem. What does metallization approach pre- 
dicts for the same problem? Assuming a <C r, it is easy 
to write down all relevant energies: E met = Q* 2 /2Dr, 
which corresponds to the uniformly smeared net charge 
Q*\ E rand = E met + Me se if, which is the energy of the 
sphere with M = 47rr 2 ro randomly distributed Z-ions on 
its surface; and Ewe = E ran d + Me{n), which is the en- 
ergy of SCL of Z-ions. Here e(n) is the energy of WC 
per Z-ion given by Eq. (pi). The last term of E\yc 
is the correlation energy p8j , it is the negative (favor- 
able) contribution due to the difference between corre- 
lated and random configurations of Z-ions. As expected, 

Emet < Ev/C < E ran d- 

To address equilibrium charge inversion, let us now 
balance the chemical potential of Z-ions at the surface 
and in the bulk. Apart from translational entropy con- 
tribution negligible at high enough concentrations of Z- 
ions, we have: for the bulk part chemical potential is 
equal to the self-energy e se if] and for the surface part 
fi = dE/dM, where E is either E met , or E ran d, or Ewe- 
Therefore, depending on the approximation we want to 
use, the equilibrium condition reads 



metal : 


Sself 


random : 


Sself 


SCL : 


Sself 



Looking at the result starting from the correlation theory, 
we see that self-energies cancel and Q*$cl — L65Zer / R, 
producing Eq. (^|). By contrast, random distribution 
leads to Q^and = 0' i- e -' charge inversion is impossi- 
ble without correlations. Finally, metallization approach 
yields Q* met = Zer/2a. Comparing with result for SCL, 
metallization approach is off by a numerical factor only 
when R ~ 2a, while it substantially overestimates charge 
inversion at R a. 

We can repeat all the above arguments for adsorption 
of PE of length L < r, radius a, and charge density be- 
low OM threshold, r\ < rj Cl on a sphere with radius r and 
charge — Q, At the surface of sphere PE molecules form 
one of several possible anisotropic WC-like phases. For 
example, it can be nematic SCL with distance A between 
almost parallel neighboring PE molecules or a solenoid 
with distance A between turns, where PE molecules 
continue each other simulating a long PE which winds 
around the sphere. In any case, balancing chemical po- 
tentials leads to the equations similar to Eqs. (|l2|): 



metal 
random 
SCL 



eseif = Q* met r}L/Dr 

£seif = Q* and r)L/Dr + e self 

Sself = Q*SCLV L /Dr + e se if 



(r, 2 /DL)ln(L/A) , (13) 



where e se if = rfLD^ 1 \n(L/a). Thus, in the correla- 
tion theory, we get Q* SCL = (ri] 2 L/Ze 2 )\n(L/A) be- 
cause of cancellation of self-energies. Random distribu- 
tion of positions and orientations of PE molecules does 
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not lead to charge inversion. For metallized state, we 
have self-energy only in the bulk, which gives Q* met = 
{rrj 2 L I Z e 2 ) ln(L / a) . In complete analogy with the pre- 
vious case, we see that metallization approach generally 
overestimates charge inversion, particularly at A a, 
i.e., when few PE molecules are adsorbed. With grow- 
ing coverage the metallization approximation becomes in- 
creasingly accurate. 

Qualitatively new effect, however, becomes important 
when macroion is so strongly charged that its neutraliza- 
tion requires almost full layer of Z-ions. Excluded vol- 
ume effect of hard cores of Z-ions adds a positive term 
to the chemical potential of SCL, which is proportional 
to fcgT and diverges at the full coverage simultaneously 
with pressure. Close to the full layer this term compen- 
sates and then over-compensates the negative Coulomb 
term /iwc-, so that charge inversion disappears. In the 
language of images, this happens because a full layer is in- 
compressible and does not allow for correlation hole and 
image. For even larger macroion charge, the second layer 
starts to form, launching a new wave of charge inversion. 
In the beginning, charge inversion is small because all 
the attraction of a new Z-ion approaching the surface is 
provided by a weak interaction with an inflated image 
in the emerging second layer, where once again A 3> a. 
Continuing, we arrive at the prediction of oscillating in- 
verted charge Q* as a function of Q (see Fig. 6), where 
charge inversion vanishes every time when the top layer 
of Z-ions is full pj|. Metallization approximation fails 
to capture the physics of these oscillations. 




1 2 



Q/Q 

FIG. 6. Schematic plot of the oscillating inverted charge 
of the sphere Q* as a function of the absolute value Q of its 
bare charge. Here Qo is the charge of the full layer of Z-ions. 
Dashed line obeys Eq. (^), which is derived for Z-ions of 
vanishing radius a. 

The oscillations of charge inversion we arrived at are 
similar to the oscillations of the compressibility and 
magneto-capacitance in the quantum Hall effect, which 
are related to consecutive filling of Landau levels |3(| . In 
this sense, we deal with a classical analog of the quantum 
Hall effect. 

We would like to conclude this subsection emphasizing 
the difference between metallization approach and the 
PB one. PB approximation effectively smears Z-ions ev- 
erywhere, both at the macroion surface and in the bulk 
of solution, while the metallization approach keeps Z- 



ions of the bulk discrete. Therefore, it is not surprising 
that the metallization approach somewhat overestimates 
charge inversion, while there is no charge inversion in PB 
approximation. 



B. Counterion release 

Consider now adsorbtion of DNA, or other polymer 
charged above the OM threshold r\ c . The result of cor- 
relation theory for this case is presented in Sec. III. An- 
other, seemingly independent reason for charge inversion, 
was suggested in jfi), arguing that the entropy gain of 
counterions released by DNA drives charge inversion. Of 
course, this theory does not pretend for universality, as 
it applies to highly charged polymers only. Nevertheless, 
we want to understand the relations of counterion release 
and correlations theories. 

Let us start with a simple qualitative argument. Imag- 
ine that neutralizing DNA molecules, along with their 
OM-condensed small ions, are adsorbing on the macroion 
(possibly releasing some of the ions), and currently the 
neutralization condition is achieved. Assume further that 
DNA rods are distributed randomly, uncorrelated in both 
positions and orientations. In this case, next arriving 
DNA molecule feels no average field, so that it has no 
reason to release its counterions. The situation is com- 
pletely different if DNA molecules are correlated on the 
surface, where locally each molecule is surrounded by a 
correlation hole - positive stripe of the background charge 
(WS cell). The corresponding field, or positive potential 
of the WS cell, may cause release of counterions not only 
at the neutrality point, but even if the surface overall 
is overcharged ("0(0) < 0). In other words, correlation 
hole, or adjustment of DNA molecules to each other, or 
image charge, or correlations (all synonyms!) is a nec- 
essary condition for both counterion release and charge 
inversion. 

To understand this better, let us re-examine the fun- 
damental physical conditions governing the adsorbtion 
equilibrium, namely that chemical potentials of both 
DNA and small ions should be the same at the surface of 
macroion and in the bulk of the solution. The equations 
expressing these conditions were derived in [jl4| (Eqs. 
(57) and (60) of that work). We re-derive them below, 
as Eq. (|l^) and Eq. <^). We deal below only with the 
case of weak screening, when r s 3> A. 

Consider first chemical potential of small ions. It must 
be the same for ions in the bulk solution, ions condensed 
on the surface of dissolved DNA rods, and ions on those 
DNA rods which, in turn, are adsorbed on the macroion. 
The two corresponding equations read: 

fcBrin(JVi,, 6 /JVi) - (2er ]c /D)ln(r s /a) 

fe B Thi(iVi )SS /iVi) = -eip(0) + e{2erf/D)\n{A/2-Ka) . (14) 
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The left-hand sides of Eqs. (jlj) is obviously the entropy 
loss of monovalent ions when they move from the solu- 
tion to the surface of DNA and the right-hand side is the 
potential energy gain per ion. In the entropic part, N\, 
Ni_ s b and Ni iSS are the concentrations of counterions in 
the bulk and at the surface of the dissolved and adsorbed 
DNA, respectively. In the energetic part, 0(0) is the av- 
erage potential of the uniformly over-charged plane with 
charge density a*, while logarithmic terms correspond 
to the potential of DNA having charge density — r\ c in 
the bulk, where it is screened by monovalent ions at the 
distance r s , and DNA with charge density —77* at the 
surface, where it is screened by other DNA molecules 
at the distance A/2n. As regards 0(0), it is given by 
the Gouy-Chapman formula 0(0) ~ — (2k B T/e) ln((r s + 
A*)/A*), where A* = Dk B T / {2i:ea*), because plane as 
a whole is non-linearly screened. Concentrations N\^ ss 
and A^i iS f, are related to the corresponding net charges: 
Ni, ss /N ltSb = (vo-V*)/(Vo-Vc)- According to Eq. (0) 
77* differs only logarithmically from r\ c . Therefore, we 
can neglect hi(Ni jSS /Ni^b) and use the same notation 
Ni yS for both concentrations. Taking into acount that 
Vc = k B TD/e, we get 



r? c ln(A*/a) = 77* ln(A/27ra). 



(15) 



Already at this stage we can look at what happens if we 
neglect correlations between adsorbed DNA rods. In this 
case, they are screened by small ions only, and the corre- 
sponding screening length is A* instead of A/2n. Eq. (|l5| ) 
then gives 77* = r/ c . Thus, if adsorbed DNA molecules 
were uncorrelated, positioned and oriented randomly, 
they would not have released counterions upon adsorb- 
tion. 

Let us return to the strongly correlated distribution 
of DNA. To complete the picture, we concentrate on the 
chemical potential of DNA rods. Corresponding equilib- 
rium condition is similar to Eq. (|]): 

£77*0(0) = K| + L( v * - i lc ){k B T/e) \n{N ltS /Ni) 

-{L 77* 2 /D) \n{\*/a) + (Lr£/D) In(r./o) , (16) 

where L is DNA length. In the right hand side, the first 
term, /j c ~ [i wc = ~{L 77* 2 / D) \xi(2ir\* / A), is the cor- 
relation chemical potential of DNA at the macroion sur- 
face. This is the energy of interaction of one DNA rod 
charged to —77* with its correlation hole which is a pos- 
itive background stripe of width A. The second term is 
precisely the entropy gain due to the release of counteri- 
ons. The third and fourth terms together represent the 
increase of the DNA self-energy when it moves from the 
bulk solution where it is charged to — r/ c and screened at 
r s to the surface where it is charged to —rj* and screened 
at A*. 

Fundamental equations ([l5]), (jl^) are identical with 
Eqs. (57), (60) of Ref B3|, which are derived there in 



a more rigorous way. The solution of these equations is 
given by Eqs. (10), (pi]), which simultaneously describe 
the correlation induced charge inversion and counterion 
release. 

Let us now continue to explore what happens if 
we imagine that DNA is randomly distributed on the 
macroion surface with respect to both position and ori- 
entation of rods, i. e. if we artificially drop all correlation 
effects. As we already mentioned, Eq. (|lf]) yields 77* = r\ c 
in this approximation. Furthermore, /i c = in the ab- 
sence of correlations, so that Eq. ( |l6| ) yields "0(0) = 
and a* = 0. Thus, assuming adsorbed DNA randomly 
distributed on the macroion surface with respect to both 
positions and orientations of rods, we arrive at the con- 
clusion that there is neither counterion release nor charge 
inversion. By contrast, lateral correlations leading to par- 
allel alignment of DNA molecules, can drive charge in- 
version even without counterion release, as this happens, 
e.g., for weakly charged PE, for which 770 < ?7 C . In this 
sense, we conclude that correlations are the driving force 
of both counterion release and charge inversion. 

Returning to the correlated distribution of DNA let 
us consider also a limit of very large r s , which is not 
very realistic, but interesting from a theoretical stand- 
point. According to the Eq. (0), DNA net charge 
77* increases with increasing r s , which physically means 
higher proportion of released counterions. To address 
this regime, we should make one step back including the 
ln(77o — 77*)/ (770 — 77c) term in Eq. (|l5|). For our purposes 
here it is sufficient just to approximate 77* by 770 in the 
Eq. (|l6|). First and third term in the right-hand side 
combine then to (Li^/D) \h(A/2tto) 1 which vanishes at 
A = 2wa and is negligible around this point. Therefore, 
the averaged surface potential in this case 

^(0) = (% - Vc){k B T/e) ln(JVi,,/JVi) + {ifjD) In(r,/o), (17) 

looks like the sum of counterion release and metallization 
terms. Interestingly, both of them are driven by correla- 
tions! 

We conclude repeating that the underlying physics of 
charge inversion is always determined by correlations. 
An explicit treatment of correlations provides regular 
and universal description of this phenomenon. Apart 
from charge inversion, the correlations manifest them- 
selves also in a number of other ways, including metal- 
like properties of the macroion surface over certain range 
of length scales and including the release of counterions 
for the case of highly charged Z-ions, such as DNA. 

T. T. N. and B. I. S. are supported by NSF DMR- 
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